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ABSTRACT 

Proton orbits are calculated in the electro- 
magnetic vacuum field of a magnetic point di- 
pole rotating with its angular velocity ST 
perpendicular to its dipole moment aj£ by nu- 
merical integration of the Lorentz-Dirac equa- 
tion. Trajectories are shown and discussed for 
various initial conditions. A critical surface 
is shown seperating initial positions of pro- 
tons which finally hit the pulsar in the polar 
region from those which finally recede to in- 
finity. 

1 . Int roduction . In this paper proton orbits will be 
shown which I obtained with one of my students (Laue and 
Thielheim, 1985 a and b) by numerical integration of the 
Lorentz-Dirac-equation in the electromagnetic vacuum field 
of a magnetic point dipole rotating with its angular velo- 
city Co perpendicular to its, dipole moment M. . Parameter 
values (model 1) are M, = lo JO G cirT and co / = 2o ^ s . 

The corresponding magnetic field strength in the polar re- 
gion of a spere of radius R = lo km is about 2 . lo^ G. 
The resulting light radius is Rj = 4.775 km. The procedures 
of numerical integration are described elsewhere (Marotzke 
and Thielheim, 1985) . 


The initial velocity of protons is zero, which is a rea- 
sonable specification since it turns out that particles 
very rapidly gain ultrarelativistic energies such that 
their orbit depends very little on their initial velocities 
as long as these are sufficiently small (i.e. nonrelativis- 
tic) . The orbits under consideration, therefore, depend on- 
ly on their initial positions. In a certain set of orbits 
these initial positions are defined by a given value of 
radius (R = lo, 2.2 and l.o in units of R T 1, equidistant 
values of°latitude (© = lo°, 25°, 4o°, 55 L6 , 7o 6 and 85° 

against the vector of angular velocity ) and equidistant 
values of longitude ( (f> Q = o°, 15°, ..., 345° against the 
x-axis, which corresponds to the direction of the dipole 
moment ft for zero phase Cot - r/R = o ). These initial 
positions are shown in figure 1. u 
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2. Symmetry Properties. 
In the given field con- 
figuration the magnetic 
vector is symmetric, 
while the electric vec- 
tor is antisymmetric un- 
der point reflection. In 
view of the symmetry 
properties of the equa- 
tion of motion, there- 
fore, protons starting 
from point symmetric 
initial positions with 
opposite initial moments 
move on point symmetric 
orbits. For this reason 
one may restrict calcu- 
lations to initial posi- 
tions on one of the two 
hemispheres . 


3. Particle Drift in the Definition of initial conditions 

Equatorial Plane. In 
the field configuration 
considered here the elec- 
tric vector is perpendi- figure 1 

cular to the equatorial 
plane z = o, while the 
magnetic vector in this 
plane is parallel to 
the z-axis, i.e.(£,^) 

= o. Thus, in the equa- 
torial plane, there is 
no acceleration along 
the magnetic field lines. 

Still there is an 
drift within the plane. 

This is illustrated by 
figure 2, the arrows in- 
dicating the drift velo- 
city. Asteriks 2 marj£ re- 
gions, where E > B and 
thus protons are torn 
out of the equatorial 
plane by the dominating 
electric field. (In this 
diagram, the arrows at 

the outer edge COrreS" E«B -Drift velocity of protons originating from the equatorial plane 

pond to initial positions 
with R q = 3) . 




figure 2 
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4. Protons starting in 
the Distant Zone. Or- 
bits of protons starting 
from initial radius 

R = lo and latitude 
@ 0 = 25° (with the ra- 
dial coordinate r given 
in the logarithmic scale 
log (r/rmin) ) are shown 
in figure 3 for At = 9o 
units of time. These 
protons move practically 
in the radial direction 
forming a narrow bundle 
slightly deflected to- 
wards the electric vec- 
tor (broken line) . An 
analogoug diagram for 
0 = 7o° is shown in 

figure 4 exhibiting a 
wider bundle of orbits. 

5. Protons starting in 

the Near-Zone. The to- 
pography of orbits is 
quite different for* 
small values of initial 
radius R , as is shown 
for R =1, 0 = 7o° and 

At - 2 (with°linear 
scale) in figure 5. Here, 
trajectories originating 
from a certain range of 
<fi do not eventually 
recede but finally hit 
the pulsar surface in 
one of the two polar re- 
gions (protons origina- 
ting from corresponding 
initial positions of the 
other hemisphere would 
eventually reach the 
opposite polar region) . 

6. Critica;! Surface. 
Obviously, for each pair 
of given initial lati- 
tude © Q and longitude 

<f) Q there is a critical 

initial radius R =R C (© 

<f> Q ) , which is tne mini-' 

mum of all initial radii 

R > R„ such that tra- 
o c 
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jectories starting from 
(R , O , <P ) run to in- 
finity. Thereby a cri- 
tical surface is de- 
fined (for parameter 
values as specified in 
model 1) which is il- 
lustrated by figure 6 in 
a perspective view (The 
largest radius is about 
R c = 2). 
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